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Abstract. 



This paper is devoted to the one-dimensional heat equation and the non-linear ordinary differential 
equations associated to it. 

We consider homogeneous polynomial dynamical systems in the n-dimensional space, n — 0,1,2,.... 
For any such system our construction matches a non-linear ordinary differential equation. We describe the 
algorithm that brings the solution of such an equation to a solution of the heat equation. The classical 
fundamental solution of the heat equation corresponds to the case n — in terms of our construction. 
Solutions of the heat equation defined by the elliptic theta-function lead to the Chazy-3 equation and 
correspond to the case n = 2. 

The group SL(2,C) acts on the space of solutions of the heat equation. We show this action for each 
n ^ induces the action of SL(2, C) on the space of solutions of the corresponding ordinary differential 
equations. In the case n = 2 this leads to the well-known action of this group on the space of solutions of 
the Chazy-3 equation. An explicit description of the family of ordinary differential equations arising in our 
approach is given. 

1 Introduction. 

The classical elliptic theta-function is a solution of the heat equation. Using the ansatz 

1 (z,2nit) = e~^ z2+r Wa(z,g 2 (t),g 3 (t)) (1) 

where a(z, g 2 , gs) is the Weierstrass sigma-function with parameters 32,93, we get a polynomial dynamical 
system on the functions g^it), gz{t) and h(t). It follows that h(t) is a solution to the Chazy-3 equation and 
$2(0, 93(t) are differential polynomials of h(t). Therefore we get a family of solutions to the heat equation, 
parametrized by the initial data of the Cauchy problem for the Chazy-3 equation (see details in pQ). 

This paper is devoted to the construction and study of a sequence of families of solutions ip(z,t;ri), 
n = 0, 1, . . . to the heat equation. For each n we get a family of solutions of the heat equation parametrized 
by the parameters of a polynomial dynamical system and the initial data of the Cauchy problem of some 
ordinary differential equation. For n = 2 we get a family which includes the one described above. 

The content of the paper is the following: In section 2 we collect the necessary facts on the classical 
one-dimensional heat equation, in particular describe the action of the group SL{2, C) on the space of its 
solutions. In section 3 for the solutions of the heat equation we study the ansatz 

ip{z, t) = e -^(*)- 2 +''W$( Zj x (t)), (2) 

where the coefficients of $(z, x(t)) as a series of z are polynomials of x(t) — (x 2 (t), x 3 (t), ...,x n {t)). This 
ansatz is closed under the action of the group SL(2, C). We present a general construction that reduces 
the problem of solution of the heat equation in this ansatz to the solution of a homogeneous polynomial 
dynamical system. In section 4 we reduce homogeneous polynomial dynamical systems of the previous section 
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to ordinary differential equations of a special type. Therefore for any n we get a family of solutions to the 
heat equation of the form ([2]). Each of this solutions is defined by a finite-dimensional numerical vector. In 
section 5 we study the ordinary differential equations obtained in the previous section and their relations for 
different n. In section 6 we study the recursions for the coefficients of the function x(t)). In section 7 
we obtain rational solutions for a series of the ordinary differential equations obtained. In the last section 
the special cases n = 0, 1, 2, 3, 4 are considered. We pay special attention to differential equations with the 
Painleve property. 

We express our gratitude to R. Kont and J. C. Zambrini for the useful discussions of results of this paper. 
The work is supported by the RFFI grants 11-01-00197-a, ll-01-12067-ofi-m-2011, RF Government grant 
2010-220-01-077, ag. 11.G34.31.0005. 



2 The general properties of the heat equation. 



The heat equation 

dt 2 dz 2 [ ' 

is linear with respect to ip an d invariant with respect to shifts of arguments: 
— for solutions ipi and ip2 the function ipx + tp2 is a solution, 
- for a solution ip(z,t) the function ip(z + Zo,t + to) is a solution for any constant zo and to, and thus 
■^tp(z,t) and M^iji{z,t) are solutions. 

Therefore without loose of generality we may consider only solutions ?p(z,t) of the heat equation that 
are even or odd functions of z regular at (z, t) — (0, 0). Any other solution will be a sum of such solutions 
up to a shift of arguments. 

Consider the problem of finding solutions to ((SJ with the initial conditions "0(0, t) — ipo(t), ip'(0,t) = ipi(t). 

For a solution of the form 

OO 

tf(*,t) = X>(*)|i (4) 



equation ([3]) takes the form 



fc=0 



00 k -i 00 k 



fc=0 k=Q 

and therefore is equivalent to the system ipk+2(t) = fc = 0, 1, . . . . Thus any solution of the form ^ 

is a sum of an even solution with the initial conditions ip(0,t) = ipo(t), ip'(0,t) = and an odd solution with 
the initial conditions ip(Q,t) = 0, ip'(0,t) = tpi(t). 
Let us consider a function of the form 

= e-^V(M), (6) 

where <p(z, t) is a polynomial of z. 

Theorem 1. The function ip(z,t) of the form ^ is a solution to the heat equation ([3]) if and only if it is 
a linear combination with constant coefficients of the function 

7^ cxp (^)) (7) 

and its derivatives with respect to z. Here c is a constant. 

Proof. If 4>(z,t) is a polynomial of degree n, set (f>(z,t) — J2k=o < l ) k{t)z k , 4> n {t) ^ 0. Substituting © into ^ 
and dividing by e ~h h ( t ) z ( we get an equality between two polynomials of z. The one in the left hand side 
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has the coefficient —^h'(f> n at z n+2 , and the one at the right hand side has the coefficient ^h 2 (f> n at z n+2 . 
Thus 

h(t) = j-^— 

for some constant c. Now the left hand side has the coefficient <p' n (t) at z n , and the right hand side has the 
coefficient — (n + ^)h(t)<p n (t) at z n . Thus 4> n (t) = c " 1 for some constant c n ^ 0. A solution with this 

(t — c) 7 

(j) n (t) is ip n (t) = (— l)"c„^pr (^ ^_ c exp ( 2 (t-c) ))- ^he otner coefficients give the relations 



(fc + 2)(* + l) 



&(*) + (2fc o +1) M^)^(t) ) , fc = 0,l,2,...,n-l 



which for fixed 0^+2 define 4>k uniquely up to the addition of a solution to 4>'^{t) = — v 2 ^ 1 ) h(t)(f)k(t), that is 
ipk(t). Thus for any polynomial 0(z, t) the corresponding solution is the sum of ipk(t) with some coefficients, 
and ipk(t) themselves up to constants are derivatives of (JTJ). □ 

The Hermite polynomials Hek{x) are defined by the relation 



d X * e ^{—) =i - 1)CXP {—) Hek{x) - 

Corollary 2. The function i\){z is a solution of the form ((6]) of the equation ([3]) if and only if h(t) = u_ c \ 
for some constant c and </>(z,i) * s a linear combination with constant coefficients of the polynomials with 
respect to z 

' -.He k (-= =). 



y/t — C y/t — C 

In the focus of our interest is the construction of a sequence of solutions of the heat equation starting 
with (J7J. Each of this solutions generates by differentiation with respect to z a new sequence of solutions. 
In the case t £ R, z £ R set 



I(t) = I iP(z,t)dz. (8) 

J — oo 

Let ip[z,t) be a solution of the heat equation. Then if 

= for to < t < h, 



dtfj 
~dz~ 



then I(t) = const on the interval (to,t±). 

Corollary 3. In the notions of theorem^ the conservation law for (|5J| holds: 

I(t) = const for t £ (c, oo). 

2.1 Symmetry groups. 



The detailed description of the Lie algebra of heat equation transform group and one-parametric sub- 
groups corresponding to the basis in this Lie algebra are given in [2]. We will need the action of the group 
SL(2,C) on the space of solutions of the heat equation. 

Let M = ^ ^ G SL(2, C) and 

/ \ i / \ 1 / —cz 2 \ , ( z at + b\ 

T{M) ^ t] = v^TT xp {WTd))^ {—d>M)- (9) 
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Lemma 4. Ifip(z,t) is a solution to the heat equation (J5J) 7 then so is T(M)ip(z,t) . 

The proof is a straightforward substitution of (JSJ) into ([3]) using ad — be = 1. 
Note r(M)l is the classical solution (JT)) of the heat equation. 

Lemma 5. The group SX(2,C) acts on i/ie space of solutions of the heat equation. 

Proof. We have r(M 2 )r(Mi)l = T (MiM 2 ) 1. Thus r(M) induces a representation of SX(2, C) on the space 
of solutions of the heat equation. □ 

3 Heat equations and dynamical systems. 

Any even (odd) regular at z = nonvanishing function if){z,t) may be presented in the form 

$(z, t) = e-i h ^ z2+r ^z 2s (f>(z, t), (10) 
where s is an non-negative integer and the function 4>(z, t) in the vicinity of z = is given by the series 

2fc+<5 

^>*) = ^ + |>(*)(^T*)!' 

where 6 = (correspondingly, 5 = 1). Note that this representation is unique, thus for any function ip(z,t) 
we may speak of the corresponding functions h(t), r(t), 4>(z, t). It follows from (JS|) that if if> solves the heat 
equation then s = 0. Further we assume S = or 6 = 1. 

Consider a function ip(z,t) in the form (|T0|) with s = and a function Denote by ji(M)h(t) 

and 72 (M) exp(r(<)) the corresponding functions in the representation in the form (fTU|) of the function 
T(M)4>(z,t). 

Lemma 6. 

*<"W« = JaTif h & + wtty *C*>«*<0> = <^hm e ' (im) 

Proof. For rp(z, t) of the form (|10|) with s = we have 

I ( \ ( 1 , /at + 6\, , c \ 2 ■ ( ok + , ( z at + b 



= j^H^t) with 0(z,t) = z ff + E fc>3 0*(*)^T for = p4^^(^D" We 

see 4>(z, t) has the from (fTT|) . thus the action on e r (*) and on /i(t) induced from the action of T(M) on tp(z, t) 
in (fTOl) is as in the lemma. □ 



3.1 General construction. 



Theorem 7. TTie function tjj{z, t) of the form (|10p is a solution to the heat equation ([3]) z/ and only if s = 0, 
r' = —(5 + and the function 4>(z,t) is a solution to 

^-4> = H2<t>-hH Q ct>, (12) 
at 



where 

'i d 2 



H2<P = ( ^"^72 + uz ) 0> 



2dz 
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Proof. Substituting (fTU|) and (|lip into ([3]) and presenting the right and left parts of this equation as series 
of z, up to elements of order 2 s + 6 we get 

(2s + S)(2a + 5- l)e r{t) z 2s+s - 2 + {z 2s+s ) = 0. 

Thus s = and up to elements of order 8 + 2 we get 

((5 + ~)/i(t) + r'(t))e r ^z 5 + (z s + 2 ) = 0. 

Thus r'(t) = —(S+h)h(t), and the last statement of the theorem is obtained directly from substituting (|TU)) 
into © and dividing by e ~^ h{t)z2+r ^ . □ 

Set x — (a?2, . . . We will assume that the variables x q have the grading — 4g, and all constants 

(that is the parameters not depending on t, z nor x) have the grading if not stated otherwise. Set 

Z 2k+S 

where z is a variable of grading 2 independent of x, and &k{x) are homogeneous polynomials of degree — 4fc. 
For example, we have $>2(x) = C2X2 for some constant C2 and $3(3;) = C3X3 for some constant C3. Therefore 
$(z; x) is a homogeneous function of degree 25. Let us remark that functions homogeneous with respect to 
some degree are often called quasi-homogeneous. We will not make any difference and further on call all this 
functions homogeneous. 

The homogeneity of x) implies 

„+i g 



"Ho* = L 3>, where Lq = ^ 2kx k - — . 

oxk 



k=2 



Let us consider in C" with coordinates (x2, ■ ■ ■ , x n +l) a se t of homogeneous polynomials p q (x) with 
degp q = — 4q, q = 3, . . . , n + 2, and introduce a homogeneous polynomial dynamical system 

—^= Pk+1 (x), fc = 2,...,n + l, degr = 4. (14) 
ar 

Set 

d 

L2§ = ^2pk+l(x2, ■ ■ ■ ,X n+ l) — 3>. 

We will need the following result: 
Theorem 8. The function $(2;; x) gives a solution to the equation 

rZ 2 $ = L 2 $ (15) 

if and only if 

S 2 = -4(1 + 2%, $3 = 2-^-<i> 2 , $ fc+1 = 2 j-<fr fe + (2fc + * ~ ^ + ^ *2*k-i, k>2. (16) 

ar ar 2(1 + 26) 

The proof follows from direct substitution of (|T3"]) into (|15[) . 

Corollary 9. Under the conditions of this theorem for a constant c such that &2(x) = CX2 we have u = 

~ 4(1+25) X2m 

Recall a function fix) is called non-degenerate if for a vector-function 4>{x) the equality (4>(x), grad f(x)) 
implies <f>(x) = 0. 
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Theorem 10. Let 3>(z; x) be a non- degenerate function of the form (|13p with coefficients satisfying (1161) . 
The function 

tt>(z,t) = $(z]X 2 (t),...,x n+1 (t)) (17) 

is a solution of the equation (TT21 if and only if the vector-function x(t) = (x 2 (t), . . . ,x n +i(t)) satisfies the 
dynamical system 

— x k (t) = p k+1 (x(t)) - 2kh(t)x k (t), k = 2,...,n + l. 
Proof. With the substitution (fT7| equation (fT2]) becomes 

fdx k , . „, , \ a$ n 

-Pk+i\?2, ■ ■ -,x n+1 ) + 2khx k j ^— = 0. 

The condition that $ is non-degenerate finishes the proof. □ 

Let us summarize the results of this section in a single theorem. Let p k {x), k — 3, ...,n + 2, be 
homogeneous polynomials of x = (x2, ■ ■ ■ , x n+ i), dega; 9 = — Aq, degpk = — 4fc, c - some constant, x) - 
a homogeneous function defined by (|13l) . 

Theorem 11 (Key theorem). Among the following conditions each two imply the third one: 

1 ) The function 

V>(z, t) = e -^(t)* 2 +rW$( z; x (t)) (18) 

solves the heat equation. 

2) The coefficients of $(z; x) are defined by the recursion 

* nX^ , ^ d * (2q + 5-3)(2q + S-2)^ , A ^ n 

with the initial conditions <& 2 (x) = cx 2: $3 = 2cps(x). 

3) The set of functions (r(t),h(t),x(t)) solves the homogeneous polynomial dynamical system 

jf = -(6 + l)h, |/ l = _^__£_ a:2) ± Xk=Pk+l{x )-2khx k , fc = 2,...,n + l. (19) 

The action L of the group SL(2,C) on solutions of the heat equation is described above. 
Corollary 12. Let ip(z,t) be a solution to the heat equation of the form (|18[) . Then 

r(M)VOM) = e-^ Wz2+?it) ^{z;x(t)) 
is a solution to the heat equation with hit) — r yi(M)h(t) ) e r ^ — r y 2 (M)e r ^ and x k (t) — ^t+d) 2k ■ 

4 Dynamical systems and ordinary differential equations. 

For a set of homogeneous polynomials P k{x) with degpfc = — 4k, k = 3, . . . , n + 2, the grading condition 
implies that for any k ^ n the polynomial Pk(x) does not depend on Xk+i, ■ ■ ■ , x n +±. We may further denote 
the polynomials p k (x 2l ■ ■ ■ , x n +i) not depending onx s ,..., x n+ i by p k (x 2 , . . . , x s -i)- 
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Lemma 13. The group of polynomial transforms of the form 

X 2 = c 2 x2, X k ^c k x k +q k (x 2 ,...,x k ^ 1 ), fc = 3,...,n + l, (20) 

where c k ^ are constants and q k (x 2 , ■ ■ ■ , x k ) are homogeneous polynomials, degq k = —Ak, acts on the space 
of homogeneous polynomial dynamical systems of the from (j!4[) . This action brings the system (|14p to 

— X k = P k+1 (X 2 ,...,X n+1 ), fc = 3,...,n + 2 
(It 

where 

Pk+i{X 2 , . . .,X n+1 ) = c k p k+1 (x 2 , . . .,x n+1 ) - —q k (x 2l . . .,x k _i). (21) 

dr 

Remark 14. ^q k (x 2l . . . , x k -i) is a short form for J2 q S;1k(x 2 , x k -i)p q +i{x 2 , x n+1 ). 

The proof of the lemma is a straightforward substitution. 

Lemma 15. The group of polynomial transforms (|20[) acts on dynamical systems of the form (|19p and brings 
a solution (r(t), h(t), x(t)) to a solution (r(t), h(t), X 2 (t), . . . , X n+ i(t)) of the dynamical system 

J t r = -{8 + \)K ±h = -#-—±^X* ±X k = P k+l{ X 2 ,...,X n+l )-2khX k , k = 2,...,n+l. 



Proof. Each transform (|20[) may be presented as a composition of transforms of the form 

X k = c k x k , c fc 7^ 0, fe = 2,...,n+l, (22) 

and 

s-2 

X k = x k , k 7^ s, X s = x s + o,qx 2 2 . . . x 3 s B ~ 2 with mj m = s, oq = const. (23) 

m=2 

Now ([22]) brings |19| to 

^ = -(<y+i)ft, £h = -h 2 - of ,° X 2 , lx k = c kPk+1 (^,...,^)-2khX k , k = 2,...,n+l. 
dt 2 dt 2(1 + 2o)c2 at c 2 c n +i 



and (J23J) brings (JT9J) to 

|r=-(*+i)^ ^=-^ 2 -^^)^ ^ = P, +1 (X 2 ,..,X„ +1 )-2^ fc , * = 2,... >B+ 1, 

where Pfc+i(X2, . . . ,X n +\) is defined by (|2~Tj) . The last statement is true because from (fT9|) it follows that 
F = x 2 2 ... x 3 °Zl , J2m=2 m irn = s, solves the equation j- f F = £ p k+1 (x)-^-F + 2shF. □ 

Lemma 16. If p s (x) does not depend on x s for some s, then the set (r(t),h(t),x 2 (t), ... ,x s -i(t)) solves the 
system 

-^r = ~(S+^)h, ^-h = -h 2 - ^ 2S) X21 -Jj. x k = Pk+i(%2, ■ ■ ■ ,x s -i)-2khx k , k = 2, . . . , s-1. (24) 

This solution along with the set of polynomials p k (x 2 , . . . , x s -i), k = 3, . . . , s, determines a solution of the 
heat equation which coincides with the solution (1181) described above. 

Thus, in the conditions of the lemma, the solution of the heat equation does not depend on the polynomials 

p s+ l, . . . ,p n +2- 



7 



Proof. If p s {x 2 , • ■ • , x n +i) does not depend on x s , then all p k (x 2 , ■ • ■ , x n +i), k = 3, . . . , s do not depend on 
x 8 , ... ,x n+ \ because of the grading, thus (j24j) is a subsystem of (fT9l) with only fc,(t), r(t), x%{t), ■ ■ ■ ,x s -i(t) 
as variables. Therefore <fr k are functions of (h(t), r(t), x%{t), . . . , x s _i(t)) because $2 and p k , k < s, are such 
functions, and are defined recurrently as such functions, and thus the solution (| 18[) of the heat equation 
does not depend on x s (t), . . . , x n +i(t) and is defined only by a solution to (fM)) . □ 

A polynomial dynamical system (TUT) is called reduced if it is defined by the set 

p k (x)=x k , fc = 3, ...,n + l, p n +2{x) = P„(a;), (25) 
where P ra (a;) is a homogeneous polynomial of degree — 4(n + 2). 



Lemma 17. i?ac/i solution of the form (| 18[) 0/ £/ie /iea£ equation may be obtained by our construction using 
a reduced dynamical system (|25[) and a solution to the system 



-^-r = -(5+ ■^h = -h 2 + x 2 , ^x k = x k+ i - 2khx k , k = 2,...,n, 

4: x n+i = P n (x 2 , ■ ■ -,x n +i) - 2(n + l)hx n+ i. (26) 
dt 

Proof. Let us start from a set of homogeneous polynomials p k {x 2 , ■ • ■ , x n +i), k = 3, . . . , n + 2. According 
to lemma \W\ we may assume that p k {x 2 , ■ ■ ■ , a; ra +i) depend substantially on x k for = 3, . . . , n + 1, thus in 
(|2H)) we may take A2 = — 2 (i+ 9^) x 2 an d recurrently 

X k = c k -ip k (x 2 , . . . ,x n+ i) + -^-q k -i(x 2 , . . . ,x k -i) for fc = 3, . . . , n + 1, 

dr 

c k -i and q k -i being defined by ([20|). The polynomial system will take the form we need. □ 

The system (|26p is defined by a number n and the polynomial P n {x). In this system r is defined as a 
function of /i(i) up to a constant of integration ro by r' = — (8 + and x 2 , . . . , a;„+i are defined by 

x 2 = h' + h 2 , x k =x' k _ 1 + 2(k-l)hxk-i, fc = 3, + (27) 
Substituting functions of h(t) into 

4: x n+i = P n {x 2 , ■ ■ -,x n+1 ) - 2(n + l)hx n+1 , 
dt 

we get an ordinary differential equation 

%,„+i(fc) = (28) 

of order n + 1 on h(t). This equation is homogeneous with respect to the grading deg h = —4, degi = 4. 
Summarizing the results of this section, we obtain: 

Theorem 18. Each solution of the heat equation of the form (1181) is defined by the set (n, P n , h, ro), where 
n is a natural number, 

P n is a homogeneous polynomial P n {x) of degree — 4(n + 2), 
h is a solution h(t) of the equation 2?p iii?l+ i(/i) = 0, 
ro is a constant. 



Corollary 19. Each solution of the heat equation of the form (|18[) is defined by a finite- dimensional nu- 
merical vector (n, P n , h n , ro), where P n is the vector of coefficients of the polynomial P n {x) and h n is the 
vector of initial data in the Cauchy problem for the equation 
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5 Ordinary differential equations, associated with the heat equa- 
tion. 



Let -Pn(y) be a polynomial of an (n — l)-dimensional argument, such that P n {x) for the graded vector 
x = (x 2 , ■ ■ ■ , x n+1 ) is homogeneous of degree — 4(n + 2). 

Denote by V n the linear vector space of such polynomials. 

Lemma 20. We have 

dim V n — p(n + 2) — p(n + 1) — 1, 
where p(n) is the number of partitions of n. 

Proof. The dimension of the space of homogeneous polynomials of degree — 4(n+l) of the variables X\, . . . , x n+ i 
with degXk = —Ak is p(n + 1). Thus p(n + 2) — p(n + 1) is the dimension of the space of polynomials of 
x\, . . . , x n +2 of degree — 4(n + 2) that can not be presented in the form x\p, where p is a polynomial of degree 
— 4(n +1). Taking away the polynomials p{x\, . . . , x n+2 ) = cx n + 2 we get V n . □ 

Consider the equation (|2"5)l . In the case P n = set 2?o,n+i(^) = P'n+iih). 

Theorem 21. Set 



V 1 {h) = (j t +h)h. 



Then for each n ^ 2 the formula holds 



V n {h) = {j t +2nh)V n _ l {h). (29) 



Proof. The equation 2? n +i = is obtained from 



dt 



- = -2(n + l)hx n+ i 



by the substitution (f2~T|) . For n = there are no substitutions to make and the equation is h' + h 2 =0. 
Advancing from T> n+1 to T> n+2 is a substitution of x n+2 = x' n+1 + 2(n + l)hx n+ i into 

^ = -2(n + 2)te n+2 , 

that is the action of + 2(n + 2)/i) on V n+1 (h). □ 
Corollary 22. For n > 1 we have 

V n (h) = h,™ + n{n + ljh^-^h +■■■ + 2"- 1 ■ nl h n+1 . 
Theorem 23. For the {n — 1)- dimensional vector-function 

V(h) = (2>i(fc),...,2? n _i(fc)) 
and the polynomial P n (y) € £/ie formula holds 

V P „,n+l(h) = ©„+l(/l) - Pn(V(h)). 

The proof follows from (|2l)]) using theorem |2"T1 
Set 

°° d 
d = Vj/fc+i^— • 
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Theorem 24. For P ra (x) € V n the formula holds 

j t +2{n + 2)hj V Pn . n+1 (h) = V Pn+un+2 (h), 

where 

P„+i(y)=dP n (y). 
Proof. For any monomial Cj ll ...j„£>i(/i) J1 . . .V n (hy n of P„(2?(/i))) we have 

A + 2(n + 2)/ij c ju ..., jn V l {h)^ ...V n (h) j » = 

= c ju ... Jn J2jk^i(h) jl • ■■■Dkihy*- 1 . ..V n (h) j » ( ^ + 2(k + l)h) V k (h) 



, Y^.iiPvhy . ..v k {hy k ~ x . ..v n (hy-v k+1 (h). 



Thus 



^- + 2(n + 2)h ) P.APUnu (,//>„)( D( h) ). 



□ 

For rt > 1 the polynomials P n may be nontrivial and we get families of equations with coefficients of 
this polynomials as parameters. For fixed values of this parameters we get special classes of differential 
equations. Further we will consider in detail the equations T>p n ^ n+ i(h) = 0, where n = 0,1,2,3,4, and 
discuss the resulting classes of equations. 

An ordinary differential equation is said to have the Painleve property if one can construct a single- valued 
general solution of this equation. A characteristic of an equation with the Painleve property is that its critical 
singular points do not depend of the initial data (see [3]). Equations with the Painleve property play an 
important role in modern physical problems. Such problems related to the Chazy-3 equation are described 
in 0], [S], [B]. The action of the group SL{2, C) on the space of its solutions is studied in [7J. Further we will 
show that for some parameters the equation Pp 2 ^{h) — coincides with Chazy-3 and Chazy-12 equations, 
thus has the Painleve property. 

According to corollary there is an action of SL(2, C) on the space of solutions of each the differential 
equation on h we cite as T>p n . n+ i(h) = 0. 

6 Discrete differential equations associated to the heat equation. 

In our general construction ( section 13. ip we construct homogeneous functions $(z; a?) that satisfy the 
equation (|15[) and thus provide solutions to the heat equation (see theorem ITT]) . In this section we return to 
the problem of construction of such functions. 

The homogeneous function $(z; x) (see ((El)) can be presented in the form 

^ 7 Z \\J\\+& 

$(z;x)=Z 5 + < J ) X fllJll+aV ' (30) 

where as before x = (x2, • ■ • dega^ = — 4fc, degz = 2, J = (j2, . . . , j n +i) is a multiindex, a(J) are 

constant coefficients and ||J|| = J2k=2 ^kjk- That is in the notions above $k(%) = j\\=2k a {J) xJ ■ 
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For a reduced dynamical system (p?5]) the equation (fT5l) has the explicit form 

__ + „/)* = P„ w _* + g> +1 _*. pi) 

We have u = — ^^s) X2 w ^ ieie c is a constant such that $2 = cx 2 . Note that in the previous sections we 
used c = —2(1 + 28). Let us present the homogeneous polynomials P n (x) in the form P n (x) = ^2 S p(S)x s 
with a multiindex S. We have \\S\\ = 2(n + 2). 

Example 25. For n — 1 we have J — (j 2 ) and the coefficients a(j 2 ) are defined by the recursion 

a(h) = 2KT2S) iih + 6 ~ 3)(4j2 + 6 ~ 2Hj2 ~ 1} 
with the initial condition a(0) = 1. 

Example 26. For n = 2 we have J — (j 2l J3) and the coefficients 0(^2^3) o,re defined by the recursion 

a(J)= 2(1 ^ 2(5) (||J||+^-3)(||J||+^-2) a (J-(l ; 0)) + 20 2 + l)a(J+(l,-l)) + 2(j3 + lM2,0) a (J+(-2,l)) 

wi£/i £/ie initial conditions a(0, 0) = 1 and a(j 2 ,j 3 ) =0 i/j'2 <0 or J3 < 0. 

Tia(j 2 ,j 3 ) = a(j 2 - 1, j 3 ), T 2 a(j 2 ,j 3 ) = a{j 2 + IJ3 - !)■ 
XTiws i/iis recursion is equivalent to the discrete differential equation Wa(j 2 , j 3 ) = /or i/ie operator 

W = 2{l C +2S) (\\J\\ + S - 3)(|| J|| + 5 - 2)7\ + 2(j 2 + 1)T 2 + 2(j 3 + 1>(2,0)T 1 T 2 - 1 - 1 

wzt/i t/ie initial conditions {a(0,0) = 1, a(j2, J3) = for min(j2,j3) < 0}. 
Set 

Tia(J) = a(j 2 - l,j 3 , ■ ■ ■ ,jn+i), T k a(J) = a(j 2 , . ..j k + 1, j k+1 - 1 . . . , j n+1 ), k = 2,...,n, 

and 

T s a{J) = T «+ s 3+-+«»+i-i T «3+-+ S »+i-i _ _ _ T -»+i-i (j) = Tf 1 ^" 1 . . . T~ 1 a(J - S). 

Theorem 27. TTie coefficients a( J) are defined as a solution to the discrete differential equation Wa( J) = 
for the operator 

n 

W = ° A \\J\\ + 5 - 2)(|| J\\ + 6 - 3)Tl + Y, 2 Uk + l)T k +2j2Un+i + 1)P(S)T S - 1 

^ ' k=2 S 

with the initial conditions {a(0, 0) = l,a(J) = for min(jfc) < 0}. The equation Wa(J) — expresses a(J) 
as a linear combination of a(J') with \ \J'\\ < \\J\\- 

The proof is a direct substitution of (|30|) into (|3Tj) using ||5|| = 2(n + 2). 

Corollary 28. For c ^ and p(S) ^ for any S we obtain a(J) /or cmy J. 

A Hurwitz series over commutative associative ring A is a formal power series in the form 



<p(z) = Y Vk-rr G -4 ® Q[[z]] with <p k e A. 



Corollary 29. If ^_^ 2 S) ^ ^ an< ^ P(^) G Z, £/ien ^(z; cc) is a Hurwitz series of z over the ring 1\x\. 
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7 Rational solutions. 



In this section we describe a construction of differential equations that leads to examples of rational 
solutions of differential equations £>p n , n +i = 0. 

For some function h = h(t) and a constant b consider an (n + 2) x (n + 2) -matrix 



Sn(h) 



bh 
bh' 



-1 
bh 



n! 

.^ (n+1) 



(n-1) 



/j(n-l) 






-(n+1) 
6/i 



\ 



The equation 



detSWW = 



(32) 



is an ordinary differential equation homogeneous with respect to the grading deg h = —4, deg t = 4. We have 

1 



f detSofW = h' + bh 2 
b 



and 



i det S B (/i) = /i (n+1) + (n + 2)6/i/i<") + • • • + b n+1 h n+2 



for n > 0. 



Theorem 30. For any integer n f/ie function 



n+l 



ft(t) = ft„(t) = i ^ -J— 

h ^ t - a k 



(33) 



fe=l 



zs </ie general solution to the equation l|32[). 

Proof. Consider the Newton polynomials s 9 = z fc of n + 1 variables where Zfe(t) = u^-t) • ^ e 

have ^-h^'it) — s g+1 . Using a classical result of the theory of symmetric functions we get S n (h) = (n + 
2)!(T n+ 2( 2; i: • • • 7 ^n+li 0) = 0> where a n+ 2() is the elementary symmetric function. Thus the function (|3"3"| is 
the solution to the differential equation (I3"2l of order (n + l). It defines the general solution to this equation 
because it depends of (n + 1) parameters a%, . . . , a n+ i. □ 

Corollary 31. For each positive integer n for the existence of a polynomial P n such that 

fp„, n+ i(/i) = -r det S n (h) 



it is necessary that 

b= (n + l). 

The calculations for n = 0, 1, 2, 3, 4 show that this condition is sufficient (see below). 
Problem 32. Prove that in the notions of the corollary the condition b = (n + 1) is sufficient. 

8 Special cases. 



In this section we construct explicitly in the cases n — 0, 1, 2, 3, 4 the formulas we have derived above. 
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8.1 Case n = 0. 

We have T>i(h) = h' + h 2 . In this case there are no variables Xk and Pq = 0. The system (j2l))) has the form 

dt K 2' ' di 
The general solution of this system is 

where the parameters rg € C and (a : b) € CP 1 are constant. Thus the function 

expro / az 2 \ 5 



solves the heat equation. In the case 5 = 0, ro = 0, b/a = c it coincides with the solution ([7]) of the heat 
equation. 

Corollary 33. For 6 = and b/a < t < oo the function ^jj{z,t) is the Gaussian density distribution. 

8.2 Case n = 1. 

We have 

V 2 (h) = (l+4/i)Di(/i) 

and deg£ 2 = —8, degPi = —12, thus Pi (2:2) = 0. The system (|26|) has the form 

d /c-1\j d ,2 ^ „, 

— r = -(d + ~)h, -r; h = - h + x 2, —x 2 = -Ahx 2 . 
dt 2 dt dt 

Therefore h{t) is a solution to 

h" + 6hti + Ah 3 = 0. (34) 
Remark 34. Differentiate (|34[) and pui y(i) = 2h(t) to get Chazy-4 equation: 

v'" = 

The general solution to (fM)) has the form 



y"' = -3yy" 3y' 2 - 3yV. 



We have x 2 (<) = -± - ^) . 

In the case a = b we come to the case n = 0. 
According to the key theorem, we have 

<S>(z;x 2 ) =/ + ]T$ 2g (:r 2 )— -— , (36) 
t/— 1 ^ '* 

where 

$ 2 = -2(1 + 2S)x 2 , <S> 2q = -(4g + S - 3)(4o + <S - 2)x 2 $ 2g -2, g = 2, 3 . . . . 
Note the function $(z; x 2 ) solves 

d 2 

^$(z; .t 2 ) = -x 2 z 2 $(z; x 2 ). (37) 
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Lemma 35. <5>(z;x2) = z 7(z 4 ; X2, 5), where , y(v,X2,S) is a solution to the differential equation 

Av 1 

7>) + 3T^ 7>) = A7(w) ' mth X = ~WTW) X2 

with the initial condition 7(0) = 1. 

Thus the function 7(1") is the eigenfunction of the generalized shift operator, defined by the generator 

, d 4 d 2 

d = 



dv 3 + 2(5 dv 2 ' 

8.3 Case n = 2. 

This case leads to remarkable equations with Painleve properties (see, e.g., [3]). 
We have P^ixi, £3) = c A x\, where C4 is a constant and 

W) = + 6h)V 2 (h), Vp 2 . 3 (h) = V 3 (h) - aV^h) 2 . 
dt 

The system (|26|) has the form 



It r = - {S+ 2 )k > dt h = - k+X2 > dt X2 = X3 ~ 4hX2 > dt 



Therefore h(t) is a solution to 

hi" + \2hh" - 18(/i') 2 + (24 - c A )(h' + h 2 ) 2 = 0. (38) 
Equation (f38|) is brought by the substitution = —6h(t) to 

y'" = 2yy" Z{y>) 2 + ^=^(6^ - y 2 ) 2 . (39) 

^ or ^ItF 4 = ~ fc^-36 e( l ua tion (fBTJf is Chazy-12. For C4 = 24 this equation is Chazy-3. 
For C4 = 6 the equation ([39]) takes the form 

y'" = 2yy" y 2 y' + ^y\ (40) 

The equation (|4l7)) is linear with respect to derivatives. 
We have 

i det S 2 (h) = h!" + Abhh" + Zb{h') 2 + 66 2 /i 2 h' + b 3 h 4 . 



For 6 = 3 we obtain the equation (|38|) with C4 = —3. 
Corollary 36. TTie function 



111 



i/ie general solution of the equation Chazy-12 (see (|39[) ) ura'i/i C4 = —3 (7c 2 = 4,). 



We will describe solutions of the heat equation in the case when y(t) = —6h(t) is a solution to Chazy-3. 
Such solutions may be constructed explicitly in terms of the Weierstrass sigma-function. 
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8.4 Necessary facts from the theory of elliptic functions. 

Details of the basic facts from the theory of elliptic functions see e.g. [3]. 
Consider the elliptic curve in standard Weierstrass form 

V = {(A, fx) G C 2 : fi 2 = 4A 3 - g 2 X - g 3 }. 

It is non-degenerate for g 2 ^ 27<7§ . Set 

f dX f XdX 

2uj k =(h — , 2r)k = -f , fc = 1, 2, 

Ja k M Ja k M 

where ^ and are basis holomorphic differentials and a k are basis cycles on the curve such that 

■qiLU 2 - = y- 

A plane non-degenerate algebraic curve V defines the rank 2 lattice T C C generated by 2co\ and 2uj 2 
with Im^ > o. 

The Jacobian of the curve V is the complex torus T = C/T. 
An elliptic function is a meromorphic function on C such that 

/(« + 2wi) = /(*), f(z + 2uj 2 )=f(z). 

The Weierstrass function p(z; 32,53) is the unique elliptic function with periods 2wi, 2ll>2 and poles 
only in lattice points such that 

H (p<*> - h) = °- 

It defines the uniformization of the elliptic curve in the standard Weierstrass form 

p'(zf = 4p(z) 3 - g 2 p(z) - g 3 . 
The Weierstrass cr-function is the entire odd function a(z) = a(z;g 2 ,g 3 ) such that 

(\ncr(z;g 2 ,g 3 ))" = -p(z;g 2 ,g 3 ) and lira (^^j = 1 - 

Periodic properties: 

a(z + 2uj k ) = —a(z) exp (2r] k (z + uj k )) , k = 1,2. 

Degenerate case: 



4 4 8 6 \ 1 1 2 2 

it I 2:; —a , — a = — exp —a z I sin az. 
' 3 27 / a P \6 f 



Consider the fields on C 2 



, . d d d 1 2 d 

'0=4.92^ ^653 — , « 2 =6g 3 - ■ 

032 033 dg 2 3 dg 3 

We have [J 0) l 2 ] = 2Z 2 , Z A = 12A, l 2 A = 0, (Z , la) = |A. 

Theorem 37 (The Weierstrass theorem). XTie operators 

Qo = zy z -l-lo, Q^i^ + i^-fa 
annihilate the sigma- function, that is 

Qoo-(z;g 2 ,g 3 ) = 0, Q 2 cr(z; g 2 , g 3 ) = 0. 
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Theorem 38 (See [I]). The function ip(z,t) such that 

1>(z, t) = e -iHt> 2 +rW a {Zi g2 ( t)) g3{t)) (41) 
for some functions r(t), h(t), 52W and g 3 (t) satisfies the heat equation 

|^) = ~|^M) (42) 

if and only if the functions r(t), h(t), g2{t) and g 3 (t) satisfy the homogeneous polynomial dynamical system 
in C with coordinates (r, h,g2,g 3 ), dcgh = —4, degr = 0: 

r ' = -\ h > h' = -h 2 + -^g 2 , g 2 = 6g 3 -4hg 2 , g' 3 = ±& - 6hg 3 . (43) 



Example 39. The function 



, ., , smaz 
ipyz, t) = exp — a t , a = const 



is a periodic odd function of z with the initial conditions i/j(0,t) — 0, ijj'(0,t) — exp (— ia 2 i). It is a 
classical solution of the heat equation. In this case we have 

1 2 , 1 2 4 4 8 fi 

r = -^t, h=-a, 92 = p, 33 = ^7- 

Example 40. For the classical solution 

1 / (z - af 
ip(z,t) = ip*(z,t) — tp*(—z,t), where ip* (z, t) = —= exp 



Vt ' V a 

which is decreasing when z — ► ±oo, we have 7 = — and 

a 2 -it , / 2a . a 2 4a 4 8 a 6 

k = -^^> r = hl( 7P )_ 2i' 52= 3F' g3 = ~27¥- 

For this solution we have ip(0, t) = and ip'(z, t) — exp ^— ^ ■ 

Consider the dynamical system (|43|) . It corresponds to the system generated by the field l 2 - 

d r d 1 2 f a a\ 

^Tff2 = 6g 3 , —g 3 = ~g 2 . (44) 

For x 2 — 72.92, x 3 = 553 we get the case n = 2, S = 1, C4 = 24. 
In this case the equation (|3"51) takes the form 

h!" + 12hh" - I8(h') 2 = 

and is brought by the substitution y(t) — —6h(t) to the Chazy-3 equation 

yi" = 2yy"-Z{y>) 2 . 

The system (j44j) has the solution 

52 (r) -3p(r + d;0,6 3 ), g 3 (r) = ip'(r + d; 0, 63), (45) 

where 63 = ^^(O) 3 — 4g 3 (0) 2 , and d is the solution of the compatible system p(d;0,6 3 ) = ig 2 (0), 
p'(d;0,6 3 ) = 2 5 3(0). 

In the terms of our construction, x 2 = \p{t + d; 0, b 3 ) for n = 2, C4 = 24. Further we will show that for 
n = 3 solutions to the heat equation with x 2 — \p{j + d; 6 2 , ^3) appear. 

The corresponding function $(z; a;) can be considered as a series of the parameter 62 and at &2 = it 
coincides with the Weierstrass sigma- function (see case n = 2). 
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8.5 Case n — 3. 



We have P 3 (x 2 , 2:3, £4) = C5X2X3, where C5 is a constant and 

V P:iA {h) = V 4 (h) - c b V 1 (h)V 2 {h). 

The system (J^SJ) has the form 

d 1., (i, , 2 d d d 

— r = -^ + -)h, — h = -h+x 2 , —x 2 = x 3 -Ahx 2 , —x 3 = x 4 -6hx 3 , — x 4 = 05X3X3 - 8ftx 4 . 
dt 2 at dt at at 

Therefore h(t) is a solution to 

h"" + 2Qhti" - 2Atih" + 96h 2 h" - UAh{h'f + (48 - c 5 )(h' + h 2 ){h" + Qhti + Ah 3 ) = 0. (46) 
Lemma 41. For C5 = 2c4 the formula holds 

zW0= (^ + 8/l ) ^2,3 

Thus for C5 = 48 the equation (j46|) is obtained from the Chazy-3 equation by the substitution y(t) — 
—6h(t) and the action of the operator (-£ + 8/1). 

Lemma 42. For C5 = 24 and y(t) = —2h(t) the equation (1461) takes the form 

y"" - lOyy 1 " + 30y 2 y" - 30y 3 y' + 6y 5 = 0. (47) 
The equation l|47p is linear with respect to derivatives. 
We have 

i det S 3 (h) = h"" + hbhti" + lObtih" + 10b 2 h 2 h" + lbb 2 h{h') 2 + 10b 3 h 3 ti + b^h 5 . 



For b — 4 we obtain the equation (j46|) with C5 = —16. 
Corollary 43. The function 



1/1 1 1 1 



4 yt — ai t — a 2 t ~ a 3 t — 04 
is the general solution of the equation (|46[) wzt/i C5 = —16. 

Consider the functions (|45|) (case n = 2) as the limit of functions 

g 2 {r) = 3p{r + d; 6 2 ,6 3 ), 93(t) = ^p'(r + d; 6 2 ,6 3 ) (48) 

at 62 — ¥ 0. The functions (14"51) along with (74 = give the general solution to the system 

d d 1 o ^ ^ 

^£2 = 633, —g 3 = -,o 2 + 2g 4 , ^-.94 = 0. 

For x 2 = j^g 2 , x 3 = |# 3 , x 4 = + .94 this system 

d d d 

—x 2 =x 3 , — 2:3 = 2:4, —2:4 = -962:22:3, 
dr ar ar 

corresponds to the case n — 3, C5 = 48. 

Therefore the solution ip(z, t) of the heat equation corresponding to this case has a parameter 6 2 , and for 
6 2 — > it tends to the solution corresponding to the case n — 2, C4 = 24. 

Thus the function e5 h (*> 2 - r W^(z;t) at t = is a deformation of the Weierstrass sigma-function 
<j(z; g 2 , g 3 ) with the deformation parameter 62. 
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8.6 Case n = 4. 



■PrAh) [j t + 10h)D,(h). 



We have P 4 (a;2, X3, £4, x$) — cq i2 x 2 + cq^x^ + £6,4X2X4, where Ce,2, C(,,z and C6,4 are constants, 
V Pi>5 (h) = V 5 (h) - c^V^hf - c 6<3 V 2 (h) 2 - <% A Vx{h)V3(h). 
The system (j2"6")l has the form 

— r = — (0 H — )n, — h——h + x 2 , —x 2 =x?—4:hx 2 , —x 3 — XA—bhxs, 
dt 2 dt ' dt A 1 dt 6 6 ' 

d d 32 

— X 4 = 2:5 - 8flX 4 , —X 5 = C 6 ^ 2 X 2 + C6,3^3 + C 6A X 2 X4 - Whx 5 . 

Lemma 44. For C6,2 = 0, C6,3 = C5 and C6,4 = C5 £/ie formula holds 

VptAh) = (j t + 10h ) v p*A h )- 

Thus for C6,2 = 0, C6,3 = 48 and C6.4 = 48 the equation T>p A ^(h) = is obtained from the Chazy-3 
equation by the substitution y(t) = —6h(t) and the action of the operator 

4 + Wh] ( 4 + 8/1 



Lemma 45. for C6.2 = — 120 ; C6.3 = 24, c 6 4 = 44 and y(t) = —10h(t) the equation T>p i5 (h) = takes the 
form 

y'"" - iyy"" + 3y 2 y>" - ^y" + AyV - J^ 6 = 0. (49) 
The equation (|49p is linear with respect to derivatives. 
We have 

det Si(h) = h'"" + 6bhh"" + 15b(ti + bh 2 )ti" + 10b(h" + Gbhti + 2b 2 h 3 )h"+ 

+ 156 2 ({h'f + Zbh 2 ti + b 2 h 4 ) h! + b 5 h 6 . 

For b — 5 we obtain T>p A ^(h) with c§ y2 = — 45, C6,3 = —26 and C6,4 = —31. 
Corollary 46. TTie function 

1/1 1 1 1 1 



5 \t — a\ t — a 2 t — 03 t — a 4 t — 05 
is tte general solution of the equation T>p i ^{h) = with Cq^ 2 = —45, cg,3 = —26 and C6.4 = —31. 

9 Addendum. 

In the main part of the paper we have examined the ansatz (|18[) for solutions of the heat equation ([3]). 
Remark that ansatz (| 18[) is a special case of 

if>(z,t) = e r{t) ^{z-x{t)) (50) 
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where 

z 2k+8 

*(*;*) = *' + E**(*)pj^j, (51) 

deg z = 2, a; = (xi, x 2 , . . . x n +i), deg x 9 = — 4q, and ^fc(x) are homogeneous polynomials of degree —4k. 
Using the approach developed in sections 3 and 4, we come to the following results: 

Theorem 47. Among the following conditions each two imply the third one: 

1) The function tp(z,t) in the form f|50f) solves the heat equation. 

2) The coefficients of i S(z;x) are defined by the recursion 

n+l g 

^ k+1 {x)=2^p j+1 {x) — ^ k {x) + ^ 1 {x)^ k {x) 

and the condition 



dxj 

J= i j 



r'(t) = ^ 1 (x(t)). 



3) The functions x k (t) satisfy the homogeneous polynomial dynamical system 

d_ 

dt 



^-x k =pk+i(x), k = 1,2, ...,n+l. (52) 



(t-a 3 ) 2 



Example 48. Torn = 2 letx x {t) = \ + ^ + . Then for x 2 (t) = -± + + 

and s 3 (t) = | ( IF ^7)3 + p=^r + ™e We 

-^■xi=x 2 , ^x 2 = X3, -^-x 3 = -3(4xix 3 + 3x 2 + 18x 2 Xi + 9x^). 

denote C^2x 2 ^ + 2x3-^ - 6(4xix 3 + 3x 2 2 + 18x 2 xj + 9xf)^ + *i(x). 

Theorem^] implies that for the function ^>(z;x) of the form (|5ip where ^i(x) = —\x\ and 'l'fc(x) = 
C 1 ^\(x) the function 

y {t - ai)(t - a 2 ){t - a,3) 

satisfies the heat equation. 

Using the methods of section 4, we bring the system (|52[) to the form 

— Xfc = Xfc+i, fc=l,2, ...,n, -rx n+ i = K n (x). 
dt dt 

Theorem 49. Each solution of the heat equation of the form (1501) is defined by the set (n,K n ,xi,r ), where 

n is a natural number, K n is a homogeneous polynomial K n {x) of degree — 4(n + 2), xi is a solution X\{t) 

of the equation x^j H+1 ^ = _R"„(xi, x' 1; . . . , xj[ ) and rp is a constant. 



Corollary 50. Each solution of the heat equation of the form (|50[) is defined by a finite- dimensional nu- 
merical vector (n, K n , C n , r ), where K n is the vector of coefficients of the polynomial K n (x) and C n is the 
vector of initial data in the Cauchy problem for the equation X]™ = K n (xi, x' 1; . . . , x^). 

The linear space of homogeneous polynomials K n (x) of degree — 4{n + 2) is p(n + 2) — 1 dimensional, 
where p(n) is the number of partitions of n (compare to lemma [2"0")l . 

Corollary 51. For each n the space of solutions of the heat equation of the form (|50p is p(n + 2) + n + 1- 
dimensional and has the p(n + 2) — p(n + 1) + n + 1-dimensional subspace of solutions of the form (|18[) . 



Therefore the introduction of the solution of the heat equation ansatz described in the main part of the 
work picks out a special class of non-linear differential equations of the form x^ 1 = K n {x\,x' Xl . . . , x x ) 
(see section [5]). 
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